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1 Introduction 



A few years ago a class of 1-dimensional "filter" Cellular Automata (CA) [P,[2j 
associated with the Schroedinger discrete spectral problem (SDSP) was intro- 
duced [3]. These CA are endowed with a number of interesting properties, 
namely: 

• they possess an infinite (countable) set of constants of motion 

• they are time reversible 

• they exhibit an interesting dynamic (solitons) 

However the authors of [3] were not able to fully develop the usual integration 
scheme (solving the direct and inverse spectral problem) for these CA; this is 
due to two technical reasons: 

1. these CA are "filter" CA 2 : i.e. their evolved state at a given time t at 
a given lattice position n depends not only on the state of these CA at 
time t — 1 but also on the "evolved" state of the CA themselves at the 
positions j < n; thus, although these CA are computable, their evolution 
is not given in an explicit form; 

2. it was impossible to find a Jost solution of the related SDSP at -co, due 
to the difficulty of projecting rational expressions on a finite field. 

In this paper we introduce three new CA associated to SDSP, and therefore 
endowed with an infinite (countable) set of constants of motion, which (sepa- 
rately) overcome these two difficulties: indeed the first one (CA1) is a standard 
CA (not a " filter" one) and the second and the third ones (CA2a, CA2b) are di- 
rect projections of the Backlund Transformation (BT) of the Kac-van Moerbeke 
Lattice 4 on a finite field. This very natural way to obtain a CA associated 
to a discrete spectral problem, namely to project the related BT on a finite 
field, is usually not exploited, due to the presence of denominators in the BT 
itself: in this paper we show how this difficulty can be overcome, by selecting 
properly the modulus of the finite field and the parameterization of the fields 
to be projected. We believe that these results are important as a necessary (if 
preliminary) step in order to develop a fully integration scheme in finite fields 
(work is in progress in this direction). 

2 Spectral problem, associated CA and their 
constants of motion 

First we show how to construct CA from the Schroedinger discrete spectral 
problem, namely: 

ip(n — 1, z) + u(n)ip(n + 1, z) = (z + z _1 )-0(n, z); (la) 
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where 



n G Z; z G C; 



(lb) 



u(n) — > 1 . (lc) 

| n >oc 

Assuming that the eigenfunction ip as well as the field u depend also para- 
metrically on the the discrete time variable t (t G Z), we look for a compatible 
evolution of ip in the form: 

M 

4>(n, z,t+l) = ^2 U {j) (n, t) ip(n + 2j, z, t) . (2) 
0=0 



Even if the results exhibited in the following could be easily extended to the 
general case, for the sake of simplicity we will restrict considerations only to the 
case M = 1, — 1, = U ; in this simple case compatibility between JIJ 
and (0) yields: 

u(n, t+\) + U{n - 1, t) = u(n, t) + U(n, t) , (3a) 



u(n, t + l)U(n + l,t) = u(n + 2, t)U(n, t) . (3b) 
The ordinary solution of © obtains solving f.i. eq. (|3al ) in terms of U and 



then inserting the obtained expression of U in eq. (|3b|) . or (equivalently, even 
if the explicit solutions looks quite different), solving first eq. (|3b(l and then 
eq. (|3al) . These two solutions are the simplest Backlund Transformations of the 
Kac-van Moerbeke lattice 0] : 



n-l 



u(n, t + l)= u(n, t) + ( ^ n _ | ^ - 1 )[ b + E (u(j,t + l)-u(j,t) 



(4a) 



, / u(n+l,t) \ tt m(? + 1, t) 

u(n ,t+i)= U (n,t) +t n u { -ij+iy ^ 

where 6 is an arbitrary constant (the BT parameter) . 

In order to view eqs. as defining a CA, one should project them on a finite 
ring, say Z m = Z/m = {0, 1, .., m — 1}, thus the problem to project possibly null 
denominators arises. In -2] a strategy to overcome this difficulty was developed: 
eqs. © were considered as " modulo congruences" an in particular ea. (|3li|) 
was solved in terms of U in the finite ring Z m equating each member of the 
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equation itself to zero. Such peculiar solution reads (here and in the following 
the symbol = denotes a m-modulo congruence): 

U(n, t) = hS (u(n - 1, t + 1)) 6 (u(n + 2, <)) , (5) 

where h , although assumed to be a constant in the following, could be in 
principle an arbitrary function of the field u, while the modulo_delta function is 
defined as: 

5(x) =0 iff x ^ mod m, x £ Z m . (6) 
It admits the following simple representation: 

m—l 

6(x) = dH(x + k), (7) 
fe=i 

where d is a normalization constant; note that, if the modulus m is a prime 
number, then 

d = m - 1 => 5{0) = 1. (8) 
Inserting (|5j into (|3a|l an interesting "filter" CA obtains (CAO): 

u(n, t + 1) = u(n, t) + h 5 (u{n -l,t+l))S (u(n + 2, t)) + 

+h{m -1)6 (u(n - 2,t+l))S (w(n + 1, t)) , (9) 

whose properties were investigated in [^j . However the " filter" and thus rather 
implicit nature of CAO is disturbing in order to fully investigate the mathemat- 
ical properties of the CA itself. 

As a first result of this paper we show that, within the same strategy adopted 
in there is an alternative solution of eqs. Q that allows to construct a new 
CA (CA1) which is "standard", i.e. a non "filter" one. Indeed inserting l|3a|l 
into l)3b|l . one gets: 

(u(n, t) + U{n, t) - U{n - 1, t)) U{n + 1, t) = u(n + 2, t)U(n, t). (10) 

It is straightforward to see that the above equation admits, in the finite ring 
Z m , the following solution: 

U(n, t) = h u(n, t)u(n + 1, t)6 {u(n - 1, t)) S (u(n + 2, t)) , (11) 

where S(x) is given by (|f>!7|> and h is again, in principle, an arbitrary func- 
tion of the field (however, for the sake of simplicity, we will assume it to be a 
constant). Inserting (|ll|l in (|3a|l , the following "non filter" CA obtains (CA1): 

u(n, t + 1) = u(n, t) + h u(n, t)u(n + l,t)6 {u(n - 1, t)) 8 (u(n + 2, tj) + 

+h (m - l)u(n - 1, t)u(n, t)S (u(n - 2, t)) S (u(n + 1, t)) . (12) 
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Note that CA1 loses the symmetry 

u(n + j,t)<^u(n-j,t + l), (13) 

which implied the "time reversibility of CAO (see 0); however CA1 is parity- 
invariant and moreover it possess a countable set of constants of motion (see 
below). For the sake of completeness, let us exhibit the explicit form of CA1 in 
the simplest non trivial case (h = 1, m = 2): 

u(n, t + 1) = u(n, t) + u(n, t)u{n + 1, t) (1 + u(n - 1, t)) (1 + u(n + 2, t)) + 

+u(n - l,t)u(n, t) (1 + u(n - 2, t)) (1 + u(n + 1, i)) . (14) 

Now, as a second result of this paper, we show how to project directly the BT 
on a finite field Z m obtaining two new "filter" automata (CA2a, CA2b). 
In our opinion this direct projection should preserve all the properties of the 
integrability scheme (IIIL'1 . To do so, we need the following Theorem (it should 
be well known to the experts in number- theory, however, for the sake of com- 
pleteness and for the convenience of the reader, we give also a neat proof of the 
theorem itself). 

Theorem 1 Consider the finite ring Z m = Z/m = {0, 1, .., m — 1} with: 

m = p k , k S N, p — prime. (15) 
Consider also the set Y C Z m of all the elements: 

y s ™s- p , s = 0, 1, ... (16) 

Then: 

1) Y is a (sub)ring (mod m); 

2) its complementary X = Z m —Y is a group; 

3) (x + y ■ z) G X ', for any x G X, y £ Y, z G Z m . 

Proof. 1) It is evident that Y is closed with respect to addition and multipli- 
cation, see <|Pol) . 

2) Obviously 1 G X (see again Moreover we know from elementary 

number theory that the modulo congruence a ■ z = c admits exactly (a, to) 
different solutions z G Z m iff (a, m) divides c (here and in the following (a, b) 
denotes the h.f.c. of a and b). Since obviously (x,m) = 1 for any x G X, it 
follows that any element of A has a unique inverse in Z m . On the contrary, 
no element of Y can have an inverse ((y,m) — (s-p,p fe ) ^ 1): thus, since 
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XflF = 0, X UY — Z m , all elements of Z m having an inverse must belong to 
X . It is also obvious that X is closed (with respect to multiplication). 

3) Since it is evident that (x + y ■ z) cannot belong to Y, then it must belong 
to X (again from X C\Y = 0, XUY = Z m ) ■ 

Let us now parametrize the field u in the following way: 

u(n,t) = ci + c 2 q(n,i), (17) 

where ci, C2 are (up to now) arbitrary constants. 

As a consequence of the above Theorem, we can project directly the BTs 
P| in the finite ring Z m , m = p k , k G N, jj = prime, provided that we choose 
ci G X, c 2 G y . However, due to (|Tc|) . we have to assume: 

d = 1, q(n,t) — > . (18) 

n | > oo 

In terms of the new field q, we obtain from l|4a(l the new cellular automaton 
CA2a: 

g(n, i + 1) = q{n, t)+ 



and, from l|4b|l . the new cellular automaton CA2b 

q(n,t + 1) = q(n,t)+ 



q(n + l,t)-q(n-l,t+l) \ yr 1 + c 2 g(j + l,t) 

l + c 2 g(n-l,t+l) / W l + c 2 g(i-l,t+l)' 1 j 

Note that these two new cellular automata are again "filter" CA, however, if the 
initial datum is given on a compact support, they are computable. Investigation 
of the dynamics of the CA introduced in this paper is postponed to a subsequent 
paper. 

Since CAl, CA2a, CA2b (as well as CAO) are related to the same linear 
problem (11121) . they share the same infinite (countable) set of constants of mo- 
tion. This set was derived for CAO in ref. here we report just the first 
conserved quantities (in terms of the above parametrization): 

+ oo 

S — — 00 

K 2 = J2 9( s '°) E 9(r,0), (22) 



r=s+2 



G 



+ oo / +oo +00 +oc 

k 3 = E c 2 E «( r >°) e - E ^ r '°) 

s— — oo \ r=s+2 I— r+2 r— s+3 



(23) 
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